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Energy storage and release in dielectric materials can be described on the basis of the charge 
trapping mechanism. Most phenomenological aspects have been recently rationalized in terms of the 
space charge model Dynamical aspects are studied here by performing Molecular Dynamics 

simulations. We show that an excess electron introduced into the sapphire lattice (Q-AI2O3) can be 
trapped only at a limited number of sites. The energy gained by allowing the electron to localize in 
these sites is of the order of 4-5 eV, in good agreement with the results of the space charge model. 
Displacements of the neighboring ions due to the implanted charge are shown to be localized in a 
small region of about 5 A. Detrapping is observed at 250 K. The ionic displacements turn out to 
play an important role in modifying the potential landscape by lowering, in a dynamical way, the 
barriers that cause localization at low temperature. 
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INTRODUCTION 



SHELL MODEL AND DETAILS OF THE 
SIMULATION 



Recent theoretical and experimental efforts have fo- 
cused on the mechanism of electron trapping in dielec- 
tric materials. The space charge model based on 
the concept of polaron has been developed to explain 
trapping and its consequences from an energetic point of 
view. The main underlying idea is that electrons trapped 
into the lattice build up a space charge simultaneously 
generating a mechanical deformation (polaron). Energy 
is stored both in electrostatic (coupling) and in mechan- 
ical (internal) forms Beyond a critical concentration 
of trapped electrons, dielectric breakdown or flashover 
occurs. The mechanical part of the energy is then re- 
leased in the form of a thermal shock wave that produces 
a local heating large enough to create a plasma, eventu- 
ally leading to the rupture of the material. One of the 
main consequences is that dielectric breakdown would 
not arise because of an external driving force (e.g. a high 
voltage difference between the electrodes), but because 
of the internal stress caused by the accumulation of the 
trapped charges. Thorough investigations have shown 
that electrons can be trapped in a lattice as soon as a lo- 
cal variation of the dielectric susceptibility occurs. This 
variation can be either due to defects (vacancies, impuri- 
ties, interfaces, etc..) or to crystallographic dissymetry. 
Once the charge is trapped the lattice is distorted and 
polarized. The energy involved in this process has been 
evaluated within some realistic approximation, and it has 
been shown that an energy of at least 5 eV is stored upon 
electron trapping. 

On the other hand, trapping has been studied exper- 
imentally using the mirror method H,^. It has been 
shown that trapping can occur in pure sapphire (a- 
AI2O3 single-crystals) only for temperatures lower than 
250 K, but the trapping/detrapping dynamics in this 
kind of materials is not yet fully understood. In partic- 
ular, the dynamic behavior in the range of pico-second 
to femto-second cannot be observed in real experiments, 
while it is crucial to elucidate the charge trapping mech- 
anisms. In this work we apply Molecular Dynamics tech- 
nique (MD), which has widely proven efficiency in re- 
producing the properties of real materials |Q , to explore 
these aspects of electron trapping in a-Al2 03. Ionic po- 
larization has been described in terms of a shell model 
which is briefly outlined in the following section. We then 
show how the solvated electron can stay firmly in some 
specific lattice sites at low temperature producing a local 
polarization and how, despite the fact that the potential 
turns out to be quite deep, the electron can self-detrap 
at moderate temperatures. 



The potential used for corundum was based on the as- 
sumption of a fully ionic model Ionic polarization 
was introduced through a modified shell model. In the 
shell model introduced by Dick and Overhauser ^ , each 
ion is represented by a massive core and a massless shell 
which simulates the valence electrons. The total ionic 
charge divides amongst core and shell, the core-shell in- 
teractions being described by harmonic springs. In this 
work we use a modified shell model due to Catlow and 
Stoncham ||^. 

The assumption of massless shells is a way to intro- 
duce the adiabatic (Born-Oppcnhcimer) approximation 
in the context of the shell model. It means that shells 
have to adjust instantaneously to the present configura- 
tion of the cores. This implies that a full relaxation of 
the shells is needed at each step of the MD simulation. 
In practice, full relaxation is very difficult and costly to 
achieve, but if a less strict convergence is required, a sys- 
tematic error is accumulated during the simulation that 
unphysically damps the motion of the cores The 
slight inconsistency between force and energy translates 
into a continuous and systematic energy loss. To avoid 
this problem we adopted a Car-Parrinello-like strategy 
by assigning a small, but finite mass to the shells, and by 
treating them as dynamical variables evolving according 
to their own (fictitious) equations of motion This 
approach proved to be very efficent in the context of ab 
initio MD, and the reason is that a second order dynam- 
ics for the shells gives rise to oscillating fictitious forces 
on the cores that average out during the slow dynamics 
of the cores. We have chosen a mass of 10% the proton 
mass, so ensuring energy conservation to high accuracy. 
The Coulombic contribution to the interionic potential 
is calculated by the Ewald summation technique p^ , p^ . 
Non-Coulombic, short-range interactions between shells 
are described by the Buckingham potential |Q, i.e. 

V,.A,^(^^-^ (1, 

The sets of parameters Ay, Cy, and p for each pair 
of ionic species are displayed in Table I, and the values 
of the spring constant and the partial charges for cores 
and shells in Table II. All values have been taken from 
Rcf. ||l^, except for the Al^+ core-shell spring constant 
which was adjusted to fit the dielectric constant of a- 
AI2O3 computed from the simulation to the experimental 
value (cr = 9.8). 

The simulation has been performed on a system of 
120 ions at constant number, volume and temperature 
(canonical ensemble). We have placed the particles into 
an orthorhombic simulation box of experimental lattice 
constants, and replicated using periodic boundary condi- 
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tions. Verlet algorithm was used to integrate the equa- 
tions of motion. A time step of 0.1 fs had to be used 
in order to integrate properly the equations of motion 
for the shells, due to their small mass. Simulations were 
started with all ions at their equilibrium positions, and 
velocities taken at random from a Maxwell-Boltzmann 
distribution at the desired temperature. Before simulat- 
ing the charge trapping we have verified the accuracy 
of the model on reproducing the physical properties of 
pure Q;-Al2 03, namely the lattice energy and the mean 
square displacement. We have found a value of -160 eV at 
T — 300 K, compared to -160.4 eV from experiment p^ ] 
and -160.24 eV from Catlow and al. |15|. The mean 
square displacement also compared reasonably with ex- 
perimental measurements at 300 K and 2170 K [|l^, as 
shown in Figs. 1 and 2. These results obtained with 
our modified shell model gave us confidence to continue 
the study of charge trapping by introducing an excess 
electron into the corundum a-Al2 03 lattice. 



TRAPPING SIMULATION 

The lattice of alumina can be described as a succession 
of anionic and cationic planes with hexagonal in-plane or- 
dering. The ions form an hexagonal compact stack- 
ing, slightly distorted in order to make room for the larger 
Al^+ ions. All sites of the crystallographic arrangement 



purely classical particle. The interaction potential be- 
tween the electron and was modelled by a pure 
Coulomb repulsion at all distances. 



are non-equivalent |18|. We first run the simulation for a 
perfect crystal of 120 ions for 10000 time steps (equal to 
1 ps) to ensure thermal equilibration. To obtain a hint 
on possible trapping sites we have computed the (unre- 
laxed) potential energy surface at the (001) plane that 
contains the ions. 



Potential energy surface 

The potential of the uncharged lattice at a generic 
point (x,y) in the (001) plane is given by the Coulomb 
interaction between the ions in the crystal and a test par- 
ticle of charge unity placed at that point. Since the test 
particle is a hypothetical object, the potential has to be 
computed with all ions in the equilibrium positions of the 
uncharged lattice. The resulting surface is shown in Fig. 
3. It can be osverved the existence of local minima where 
the electron can be trapped. We have chosen the position 
of one of these potential wells as the initial position for 
the excess electron in our simulation. 



Model for the ion-electron interaction 

In this first approach to the problem of electron dy- 
namics in dielectrics we have considered the electron as a 



*e02- {r) = 



for core and shell 



(2) 



The e~-Al^+ interaction was instead described in 
terms of short-range repulsion, as though the excess elec- 
tron was an oxygen shell, but of charge unity and finite 
mass (the electron mass). 



RESULTS 

We have performed simulations at several different 
temperatures up to 300 K. Fig. 4 shows the path of 
the (classical) electron at 300 K, and Fig. 5 corresponds 
to 200 K. Also the potential energy surface in the (001) 
plane is shown in the plots. In Fig. 6 we plot the mean 
square distance traveled by the electron as a function of 
time. It is clearly apparent that at 300 K the electron 
escapes from the initial trapping site, while at 200 K 
it stays close to the bottom of the potential well. The 
transition from one regime to the other has been located 
around 250 K. Below this temperature the excess elec- 
tron solvates, and above it jumps from one well to an- 
other. These results are in good agreement with experi- 
mental data obtained by the mirror method [g . 

In no case the kinetic energy of the ions is sufficient 
to overcome the barriers present, but at high tempera- 
tures the ions in the polarization cloud can respond more 
rapidly by lowering the barriers in a dynamical way. In- 
deed, the active presence of the excess electron has an 
important backreaction effect onto the neighboring ions. 
When the excess electron is trapped, it attracts the pos- 
itive Al^+ ions and repels the negative ions causing 
a lattice distortion and a dielectric polarization. 

Despite the limited size of our system we have been 
able to determine the main features and magnitude of 
the ionic dispacements. The six nearest oxygen ions 
moved significantly outwards, away from their equilib- 
rium positions, typically between 0.17 and 0.27 A. The 
displacement pattern of the Al^+ ions is more compli- 
cated. The two first neighbors are very stable, with an 
average vibration of 0.025 A. This is because these two 
ions are strongly attracted towards the tapped electron, 
and hence they behave as though their mass were renor- 
malized by the interaction. The next two AP+ neighbors 
are also largely attracted towards the trap, because the 
outwards motion of the ions leaves enough place for 
them to move in. The average fiuctuation of the AP+ 
ions is about (0.07 ± 0.01) A. The outer AP+ ions do 
not experience any significant force, as a consequence of 
the short range of the e~-AP+ interaction. According 
to this pattern of distortions, we estimate the size of the 
polaron in about 4.8 A. 
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We have computed the potential energy surface for the 
charged lattice at T = 200 K , and we show it in Fig. 7 
compared to the uncharged potential. The former was 
calculated with the ions in their displaced positions af- 
ter the introduction of the excess electron, but only the 
potential felt by the test charge was computed, in order 
to compare with the uncharged lattice. The bottom of 
the potential well turns out to be 4.7 eV lower in energy 
in the case of the deformed (charged) lattice. Therefore, 
the effective presence of the electron further stabilizes 
the trap by increasing the depth of the potential well. 
These results are consistent with experiments showing 
that the energy necessary to detrap charging electrons is 
larger than the energy gained upon trapping them 
The total energy of the lattice increases about 4-5 eV 
upon electron implantation, so that this number repre- 
sents the energy of formation of the polaron in a-Al2 03. 
This value is consistent with the estimation done using 
the space charge model Part of the energy is stored 
as mechanical energy and, when dielectric breakdown oc- 
curs, it is supposed to be released in the form of a shock 
thermal wave that eventually leads to fracture. 



CONCLUSIONS 

In this paper we have investigated the phenomenon 
of electron trapping in the corundum (a-Al2 03) lat- 
tice. Our simulation revealed that electrons can only be 
trapped at temperatures lower than Td — 250 K. This is 
in good agreement with experimental results obtained us- 
ing the mirror method. Trapping turns out to be possible 
only at lattice sites where the potential surface of the un- 
charged lattice exhibits a minimum. This is in agreement 
with the space charge model, as these sites appear to be 
polarizability defects due to crystallographic dissymetry. 
Although our simulations reproduce the main features of 
electron trapping at a semiquantitative level, further im- 
provement is needed in two directions: (a) modelling of 
the short-range interaction of the electron with the ionic 
species, and (b) quantum treatment of the electronic dy- 
namics. Work is in progress along these two lines. 

The authors are grateful to E. Smargiassi for many 
useful suggestions, in particular on the establishment of 
the shell model. 



TABLE CAPTIONS 



Interactions 


Aij (eV) 


P (A-^) 


Cij (eV A") 


Al-Al 
0-0 
Al-0 


0* 
22764* 
1460.3* 


0.1* 
0.149* 
0.29912* 


0* 
27.88* 
0* 



Table II : Values of the parameters used for the shell 
model (* values from Catlow et al.. 0). 



Parameters 


AP+ 


0^- 


spring constant (eV A^^) 
core charge (e) 
shell charge (e) 


924.88 
1.6170* 
1.3830* 


103.07* 
0.8106* 
-2.8106* 
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Table I : Values of the parameters used for the short- 
range interactions (* values from Catlow et al. [^). 



FIG. 1. Mean square displacement for the AI^^ ions at 
T=300 K and T=2170 K. The curves are scaled by a factor 
X = 0.01 and x — 0.1, respectively. 
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FIG. 2. Mean square displacement for the ions at 
T=300 K and T=2170 K. The curves are scaled by a factor 
X = 0.01 and x = 0.1, respectively. 



FIG. 3. Potential surface for the uncharged lattice. 



FIG. 4. Path of the excess electron at T = 300 K. 



FIG. 5. Path of the excess electron at T = 200 K. 



FIG. 6. Mean square displacement of the electron at 
T = 300 /s: and T = 200 K. 



FIG. 7. Comparison between the potential surfaces of the 
uncharged lattice and the charged one. 



5 



8 




time (ps) 



